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We overview some of theory and phenomenology aspects of high energy 
factorization. In the theory part we focus on basic equations of high en- 
ergy factorization i.e. BFKL, CCFM, BK. In the phenomenology part we 
focus on forward-central jets correlations at Large Hadron Collider and on 
production of charged particles in HERA. 

PACS numbers: 12.38. Cy 

1. Introduction 

Large Hadron Collider (LHC) is already operational and Quantum Chro- 
modynamics (QCD) is the basic theory which is used to set up the initial 
conditions for the collisions at the LHC but also to calculate hadronic ob- 
servables. Application of perturbative QCD relies on so called factorization 
theorems which allow to decompose given process into long distance part 
called parton density and short distance part called matrix element. Here 
we will focus on high energy factorization [lj (there exist also collinear fac- 
torization scheme but we will not discuss it here) which applies when both 
momentum scale and energy scale involved in scattering process are high 
(for recent works on relation of high energy factorization to other schems we 
refer a reader to [2] ) . The evolution equations of high energy factorization 
sum up logarithms of energy accompanied by coupling constant a"ln m s. 
Depending on the energy range and observable one may use: BFKL [3HS]> 
BK [6, 7j or CCFM |14H16| evolution equation. When the energies of the 
collision are of the order of 10 3 GeV and one considers inclusive processes 
in electron-proton Deep Inelastic Scattering as for example at HERA the 
BFKL approximation applies. CCFM equation since it depends on the 
hardness of the probe allows additionally for studies of exclusive final states. 
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Fig. 1. Factorized structure of the cross section. 

However, if one would like to account for formation of dense system like in 
nuclei-nuclei collision where partons eventually overlap the BK equation or 
some nonlinear extension of CCFM has to be considered since it apart from 
splittings of gluons allows for their recombination. 

The paper is organized as follows. In the following section 2 we intro- 
duce framework of high energy factorization, we introduce basic evolution 
equations. In section 3 we present phenomenological applications on two 
examples: production of forward-central jets at Large Hadron Collider and 
production of charged particles at HERA. 

2. High energy factorization and evolution equations in pQCD 

The high energy factorization formula while applied to jet production in 
hadron-hadron scattering reads: 

d ° = y I d2kid2h 1 \Au^f (2.i) 
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where k\ = |ki|, k 2 = | | and x\, x 2 are longitudinal momentum frac- 
tions see FigdU). The functions A a /A(xi, k 2 , /J, 2 ) and Af,m(x 2 , k 2 , /i 2 ) are 
the unintegrated distributions which are solutions of high energy factoris- 
able evolution equations. They describe distributions of transversal and 
longitudinal momenta of partons in the incoming protons A and B respec- 
tively. This factorization formula apart from summing up logarithms of 
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energy accounts also for hardness of the process. The sum is made over all 
flavors of initial and final partons. Similar formula can be written for DIS 
process. 



2.1. BFKL equation 

The simplest of high energy factorizable evolution equations is the BFKL 
equation. At leading order in lnl/x (LLlnl/x) it reads: 



2^ 



d(j)(x, k 
dlnl/x 



a, / dl 
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l 2 (j)(x,l 2 ) - k 2 ^(x,k 2 ) k 2 (j)(x,k 2 y 



(2.2) 



where we used standard notation (j)(x, k 2 ) for factorization scale independent 
unintegrated gluon density. The real emission part of the kernel describes 
radiation of gluons which are strongly ordered in longitudinal momentum 
fraction i.e. they are well separated in rapidity. The virtual part removes 
singularity when k = I and is called the Regge trajectory of the gluon. This 
equation predicts strong rise of gluon density at small x (j)(x,k 2 ) = x~ x 
and this tendency is not changed even if subleading logarithms in In 1/x are 
taken into account |12j . Recently the BFKL equation with special ansatz for 
running coupling constant and renormalization group improved kernel has 
been used to describe F2 data and very good description of proton structure 
function has been achieved H3l. 



2.2. CCFM equation 

The CCFM is an equation which sums up cascades of gluons under the 
assumption that subsequent emissions are ordered in an emission angle. It 
assumes the following form: 

A(x,k 2 ,p) = a s [ dz [ ^-j;0(p-zq) A s (p, zq) "" ' 



7rq z \ z 1 — z 

xA[- z ,k',q) . (2.3) 

The momentum variable p is defined via £ = p 2 /{x 2 l s), and k! = |A; + (1 — z)q\ 
the momentum q is the rescaled momentum of the real gluon, and is related 
to q by q = q/(l — z). Here A s is the Sudakov form factor which regularizes 
the singularity of the 1/(1 — z) pole, while A ns is the so called non-Sudakov 
form factor and it corresponds to virtual contribution in the BFKL equation 
(their detailed form is not important in this note). 

For recent theoretical works on this equation we refer reader to [S1[S] while 
for phenomenological applications we refer reader to [10(.ll lj . The CCFM 
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equation has been derived after observation of coherence effects in emission 
of gluons |16j and it combines information from BFKL and DGLAP and 
reduces to each of them in appropriate limits: BFKL in the limit when 
x — ¥ and DGLAP when x — >• 1. Since it depends on hardness of the probe 
and also on k 2 of the incoming gluon it might be used in studies of final 
states. In particular, one can study physics of production of jets in forward 
direction which we overview in section 3. 



2.3. Saturation effects: BK equation 

As it has been already remarked, if one wants to study physics at largest 
energies available at LHC one eventually has to go beyond BFKL, CCFM. 
This is because these equations were derived in an approximation of dilute 
partonic system where partons do not overlap or to put it differently do not 
recombine. Because of this those equations cannot be safely extrapolated 
towards high energies, as this is in conflict with unitarity requirements. To 
account for dense partonic systems one has to introduce a mechanism which 
allows partons to recombine and therefore to saturate [T7j. Existing data 
suggest that the phenomenon of saturation occurs in nature. The semi- 
nal example is provided by the discovery of geometrical scaling in HERA 
data [18] and more recently by geometrical scaling in production of inclusive 
jets at the LHC data [19l[20]. Also the recently observed ridge-like struc- 
ture in p-p collision at he LHC [21] has been described within approach 
including saturation [22]. There are various ways to approach the problem 
of evolution allowing for formation of dense system |23H31| . here we are 
interested in the one which can be directly formulated within high energy 
factorization approach [I]. In this approach one can formulate momentum 
space version [32J of the Balitsky-Kovchegov equation which sums up large 
part of important terms for saturation and which is a nonlinear extension 
of the BFKL equation. The equation reads: 



(2.4) 



where R is the radius of the proton. The nonlinear term being convolution 
of the triple pomeron vertex [33] with gluon density allows gluons to merge 
apart from gluons splitting. Due to the interplay between splitting and 
merging of gluons the equation above generates dynamically scale which is 
called saturation scale Q s (equation 12.41 has been also obtained by transform 
of coordinate version of BK equation to momentum space as has been done 
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Fig. 2. Gluon density obtained from CCFM with saturation to gluon density from 
CCFM as a function of k 2 for x = 1(T 5 7 x = 1CT 6 



in |34j). This scale acts effectively as a mass of gluons and therefore regulates 
bad infrared behavior of gluon density [35]. It also selects the most probable 
k of gluon to be of order of the saturation scale. It follows from the fact 
that at k = Q s the gluon density has a maximum: 

Q s = d lnk 2^(x,k 2 ) = 0. (2.5) 

In order to have flexible approach and be able to simulate the scat- 
tering process in detail one uses Monte Carlo implementation of evolution 
equations, this is the case for DGLAP, CCFM, BFKL. As BK is a nonlin- 
ear equation it is not of straightforward usage in Monte Carlo generators. 
However, one can avoid complications coming from nonlinearity by applying 
absorptive boundary conditions [36] which mimics the nonlinear term in the 
BK equation. 

2.4- Saturation effects: CCFM equation with absorptive boundary 

The basic principle of this method is that the absorptive boundary limits 
the phase space for gluons and therefore effectively acts as nonlinear term 
in the BK equation. In the original formulation it was required that the 
BFKL amplitude should be equal to unity for a certain combination of k 2 
and x. In discussed here approach the energy dependent cutoff on trans- 
verse gluon momenta has been imposed. It acts as absorptive boundary 
and slows down the rate of growth of the gluon density. In order to have 
description of exclusive processes and account for saturation effects one can 
use CCFM evolution equation together with absorptive boundary [55] . Its 
certain variation has been implemented in CASCADE Monte Carlo event 
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Fig. 3. Production of forward-central jets. 

generator [37] (for adaptation to CCFM of method developed in [3U] we refer 
the reader to [HJ§]). In the approach of [38] as the absorptive boundary the 
GBW [39] critical line has been used. The condition for saturation was pro- 
vided by the GBW saturation scale Q s = ko(xo/x) x ^ 2 i.e. density of gluons 
with momenta generated at given x with transversal momenta which satis- 
fied condition k<Q s was set to go to zero as k 2 what is in agreement with 
numerical solution of BK. This prescription gives gluon density which has 
a maximum as a function of k in agreement with results obtained from BK. 
However, one should add that this method is quite simplistic and within 
this approach one can not find the effect of saturation of saturation scale 
itself which has been found in [9]. 

3. Phenomenological applications. 

3.1. Production of central-forward jets at LHC 

Physics in the forward region at hadron colliders is traditionally dom- 
inated by soft particle production. With the advent of the LHC, forward 
physics phenomenology turns into a largely new field [40H42] involving both 
soft and hard production processes, because of the phase space opening up 
at high center-of-mass energies. Forward jet production enters the LHC 
physics program in an essential way both for QCD studies since one can 
probe dense parton systems [43] and for new particle searches, e.g. in vec- 
tor boson fusion search channels for the Higgs boson [441145] . The forward 
production of high-px particles brings jet physics into a region character- 
ized by multiple energy scales and asymmetric parton kinematics. Here we 
overview results of [TU] where the study of forward-central jet correlations 
of two jets has been done. The results of such investigations can serve to 
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Fig. 4. Transversal momentum spectra of produced jets at total collision energy 
yfs = 7TeV with requirement that p±_ > 10 GeV. We compare predictions obtained 
from CASCADE and PYTHIA running in a multiple interactions mode and no 
multiple interactions mode. Spectrum of forward jets (left); spectrum of central 
jets (right). 



estimate the size of backgrounds from QCD radiation between jets at large 
rapidity separations for Higgs boson searches in vector boson fusion chan- 
nels. High-energy factorization allows one to decompose the cross section 
for the forward-central jet production of Fig. [3] into partonic distributions 
and hard-scattering kernels, obtained via the high-energy projectors [lj from 
the amplitudes for the process p\ + pi — > P3 + Pa + 2 massless partons. 
In Fig. 2] the prediction of differential cross section is shown as ob- 
tained from CASCADE and PYTHIA. The cross sections predicted from 
both simulations at low momentum are of the similar order, however, at 
larger transverse momentum the CASCADE predicts a larger cross section 
what is clearly visible for central jets (Fig. fright). This behavior is ex- 
pected since CASCADE uses matrix elements which are calculated within 
high energy factorization scheme allowing for harder transversal momen- 
tum dependence as compared to collinear factorization. Moreover CAS- 
CADE applies CCFM parton shower utilizing angle dependent evolution 
kernel which at small x does not lead to ordering in transverse momentum, 
and thus allows for more hard radiation during evolution as compared to 
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Fig. 5. Differential cross section for transverse momentum distribution of charged 
hadrons calculated within CCFM (violet continuous line), CCFM with saturation 
(dashed blue line) and DGLAP (dotted black line) 

based on leading order DGLAP splitting functions Monte Carlo generator 
PYTHIA. The parton shower has major influence on the side where the 
small x gluon enters the hard interaction, thus the jets in the central region 
are mainly affected by the parton shower. 



3.2. Production of charged particles at HERA 

Another exclusive observable that is interesting to look at is the pt spec- 
trum of produced charged particles in DIS. Here we review application of 
the unintegrated gluon density from CCFM with introduced saturation ef- 
fects via energy dependent cut off mimicking nonlinear effects as has been 
done in [38] . In the results we see a clear difference between the approach 
which includes saturation and the one which does not include it. The de- 
scription with saturation is closer to data suggesting the need for saturation 
effects [16]. We compare our calculation with calculation based on CCFM 
(CASCADE) and on DGLAP (RapGAP) evolution equations. From the 
plots Fig. [5] we see that the CCFM with saturation describes data better 
then the other approaches. CCFM overestimates the cross-section for very 
low x data while DGLAP underestimates it. This is easy to explain, in 
CCFM one can get large contributions from larger momenta in the chain 
due to lack of ordering in k 2 while in DGLAP large k 2 in the chain is sup- 
pressed. On the other hand CCFM with saturation becomes ordered for 
small x both in k 2 and rapidity and therefore interpolates between these 
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two. 

4. Conclusions 

We have reviewed basic theoretical aspects of high energy factorization 
and gave examples of phenomenological applications to hadron-hadron col- 
lisions and lepton-proton collisions. We also stressed the uniqueness of the 
high energy factorization as a framework in which both hard processes and 
formation of dense system can be studied. 
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